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Simulation of depositions of a Lennard-Jones cluster on a 
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Kyoto 606-8502, Japan 

Depositions of amorphous Lennard-Jones clusters on a crystalline surface are numeri- 
cally investigated. From the results of the molecular dynamics simulation, we found that 
the deposited clusters exhibit a transition from multilayered adsorption to monolayered ad- 
sorption at a critical incident speed. Employing the energy conservation law, we can explain 
the behavior of the ratio of the number of atoms adsorbed on the substrate to the cluster 
size. The boundary shape of the deposited cluster depends strongly on the incident speed 
and some unstable modes grow during the spread of the deposited cluster on the substrate. 
We also discuss the wettability between different Lennard-Jones atoms. 

§1. Introduction 

A nanocluster containing 10 - 10,000 molecules exhibits intermediate properties 
between bulk materials and individual molecules. Recently, there has been growing 
interest in the physics of nanoclusters.'^"^"^'^'-^ In particular, it is important to 
investigate depositions of nanoclusters on solid surfaces for the construction of high- 
quality films used in nanoscale electronic devices and photonic devices-l^ 

The ionized cluster beam (ICB) technique was developed by Yamada et al^^^ 
The ICB technique is used to produce atomic clusters by employing adiabatic ex- 
pansion of condensed vapour through a nozzle into a high vacuum region. In the 
ICB technique, clusters are ionized by electron impact and then accelerated toward 
a substrate. Because the ICB technique controls the translational kinetic energy 
of the cluster, there have been many experimental and theoretical studies aimed at 
understanding the influence of the incident velocities of the cluster .^3 

The outcome of such a cluster impact is largely influenced by the incident veloc- 
ity, as seen from the phase diagrams in references.13''Hli If the translational kinetic 
energy per atom becomes too large, the cluster can damage the substrate, and the 
cluster can break into pieces after the impact .IIS '^^'^^i 'US 'US} S-ED However, if the 
translational kinetic energy per atom is less than 100 eV, the cluster is adsorbed 
on the surface or reflected by the surface. Awasthi et al. carried out molecular 
dynamics simulations for collisions of Lennard-Jones clusters with weakly attractive 
surfaces .l^'l^ They discovered that the cluster rebounds when the translational 
kinetic energy of the cluster is larger than the adhesion energy. Moreover, they clar- 
ified that a transition from adhesion to rebound occurs at the critical translational 
kinetic energy. Jarvi et al. carried out molecular dynamics simulations of low-energy 
deposition of individual metal clusters on a (100) surface.lSS'lSSJ They revealed that 
the heat generated by the impact partially or completely melts the deposited cluster. 
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As a consequence, the atoms in the cluster are rearranged and adjusted to the atomic 
structure of the substrate. They found the maximum size at which single clusters 
align epitaxially on the substrate. 

Recently, Kuninaka and Hayakawa carried out molecular dynamics simulations 
of two identical colliding clusters and investigated impact phenomena of nanoclusters 
subject to thermal fluctuations.'^ They found super-rebound events in which the 
restitution coefficient is larger than 1. They confirmed the validity of macroscopic 
quasi-static theory of cohesive collisions.'^ This suggests that the research of nan- 
oclusters are relevant even for the study for fine powders whose diameters are ranged 
from 100 nm to 1 ^ j^j28li,|29} They also revealed the mechanism responsible for the 
super-rebound process, the normal rebound and the merging. 

Although early numerical studies assumed that the clusters are highly crys- 
tallised, we also need to know the properties of amorphous clusters. Indeed, it is 
easy to form metastable amorphous clusters in terms of the quench process from 
high temperature hquids.ES'EB-ElS.El-IM} 

The main purpose of our paper is to understand the behavior of the deposited 
amorphous Lennard-Jones clusters on the crystalline surface at zero temperature. 
Here, we report on our molecular dynamics simulation of the depositions with the 
small incident energies per atom which lie in the so-called soft-landing regime (0 — 2 
eV). In addition, we report on the wettability between different Lennard-Jones atoms. 

From the analysis of the final configurations of the deposited clusters, we find 
the existence of a morphological phase transition from the hemi-spherical droplet to 
the monolayer film at the critical incident speed. The multilayered adsorption state 
is described by the energy conservation law. Furthermore, we find that there are 
some unstable modes of the boundary shape of the deposited cluster. 

The organization of this paper is as follows. §2 consists of two subsections. §2.1 
explains the model of our numerical simulation. We explain our setup of cluster de- 
positions in §2.2. §3 consists of five subsections. §3.1 exhibits some time evolutions 
of impact processes. §3.2 explains how the cluster size and the cluster adsorption 
parameter depend on the incident speed after the impact. §3.3 discusses the transi- 
tion from partial wetting to perfect wetting of deposited clusters. §3.4 explains the 
morphological change of the final configuration of adsorbed atoms in clusters. §3.5 
discusses the wettability between different Lennard-Jones atoms. In §4, we discuss 
our numerical results and we summarize the conclusion. 

§2. Molecular dynamics simulation 

2.1. Model 

In order to investigate the nanocluster depositions on a substrate, we perform 
a molecular dynamics simulation. Figure [T] displays one snapshot of our numerical 
simulation. Because we are interested in neutral nanoclusters and substrates, the 
electrostatic interaction between atoms is not considered. We assume that the poten- 
tial energy of interaction between two atoms can be described by the Lennard-Jones 
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Solid surface ( 1 layer, 120 x 120 atoms) 

Fig. 1. (Color online) One snapshot of our simulation of a nanocluster deposition. The incident 
cluster contains 300 atoms which are bounded by the Lennard- Jones potential. The substrate 
consists of a single layer ( 120 x 120 ) atoms on a triangular lattice. 

( LJ ) potential: 

where subscripts a and (3 specify the species of LJ atoms, Vij is the distance between 
two atoms labeled by i and j. Here, e^f^ and are respectively the strength of 
the interaction and the diameter of the repulsive core between a atom and (3 atom. 
If a 7^ /?, we adopt the cross parameters of the LJ potential by using the Lorentz- 
Berthelot rule as 

<^af3 = ^ , = ^(^aGf^- [2-2) 

We mainly investigate the case that the cluster and the substrate are constructed 
by same atoms, A. Here, we borrow the values of LJ parameters and mass of a 
typical inert gas, argon. Therefore, eAA, ctaa and the mass of an A atom are 
1.65 X 10-2^ J, 3.405A and 6.63 x lO'^^kg, respectively.'3S''33'l3S On the other hand, 
to study the influence of the interaction energy between different atoms of the cluster 
and the substrate, we introduce a C atom as the constituent of the substrate in §3.5. 
We use the values of LJ parameters and mass of carbon to the C atom, in which 
ecc, crcc and the mass of C atom mc are 3.86 x lO^^^J, 3.354A and 1.99 x 10~^^kg, 
respectively.'^ ''^''^ We should note that the interaction energy between the cluster 
and the substrate, eAC; is several times weaker than eaa- In the following, we omit 
the subscripts of the LJ parameters of the interaction between the A atoms and we 
adopt simplified notations e, a. We also adopt m as the mass of an A atom. We 
use e, a and m as the units of energy, length and mass, respectively. Thus, the unit 
time is given by r = y^ma^/e. 

We use a single layer surface which involves 120 x 120 atoms on a triangular 
lattice as the substrate with the periodic boundary condition.'^'HB We set the lattice 
constant to 2^/^cJq,q, ( a = A, C ) as the equilibrium distance between atoms. To 
avoid the destruction of the substrate, each atom of the substrate is also tethered 
to its equilibrium position by an elastic spring. In actual impacts of nanoclusters 
on substrates, the energy induced by an impact is relaxed to the internal motion of 
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the atoms of the bulk region of the substrate. To represent such a energy relaxation 
process, we simply introduce the viscous force proportional to its velocity. The 
introduction of the viscous force has another advantage to reduce the unrealistic 
boundary effects. Indeed, if we do not introduce the viscous force, the acoustic wave 
would be transmitted across the boundary. The atom of the substrate at rj satisfies 
the equation of motion 



»"^ = -E5;-^('-.)-«r.-0-Af. (2.3) 



where is a summation over the interacting pairs i and j, and r^^ is its equilibrium 
position, and ttt-q, ( a = A, C ) is the mass of an atom of the substrate. Because we 
consider the substrate as a rigid surface, we use the spring constant k = l.Ox lO^e/o"^. 
We simply adopt the value A = 1.0^/me/a for the coefficient of viscosity. 

We adopt the velocity Verlet method for numerical integration of the equation 
of motion for each atom with the time step dt = 1.0 x 10~^y^m(j^/e. To reduce 
computational costs, we introduce the cut-off length tJcut = 3.0(7 to the LJ potential, 
and we adopt the periodic boundary conditions in the horizontal xy directions and 
the free boundary condition in the vertical z direction. It should be noted that, the 
viscous force is evaluated as the value at the previous time step for the numerical 
integration of Eq. (|2-3p . 

2.2. Setup 

We make a LJ cluster by the temperature quench^^' into the metastable phase 
of LJ fluid.H3 We prepare 32, 108, 255, 300, 500 and 862 atoms in a periodic box 
and equilibrate at the temperature T = l.Oe with the number density 0.05o"~'^ in 
the gas state ( Fig. [2] (a) ). It should be noted that the unit of the temperature 
becomes e, because we set the Boltzmann constant to be unity. To equilibrate the 
gas at a specific temperature, we adopt the velocity scaling method and perform 
until r = 2000(7 Y^m/e for the relaxation to a local equilibrium state. We have 
confirmed the equilibration of the total energy in the initial relaxation process, and 
we quench the gas to T = 0.5e. After an equilibration, a weakly bounded liquid 
cluster is formed ( Fig. [2] (b) ) and is quenched to T = O.Ole to make it rigid. This 
two step quench is adopted to form one cluster from an initial gas state. Indeed, if 
we quench the system into T = O.Ole directly, many small clusters appear. After this 
equilibration process, we obtain an amorphous cluster ( Fig. [2] (c) ). We place the 
amorphous cluster at lOci above the substrate and give the cluster the translational 
velocity Mmp to make it collide against the substrate. It should be noted that the 
amorphous cluster is metastable to keep its shape within our observation time. The 
incident angle of the cluster to the substrate normal is zero. The incident speed of 
the cluster is ranged from Vimp = 0.1-y/e/m to Vimp = ^.Q^/eJm. 
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Fig. 2. (Color online) Illustration of a LJ cluster formation consisting of 300 atoms formation by 
temperature quench method, (a) The initial configuration of atoms in the gas phase at T = l.Oe. 
(b) A liquid cluster obtained from the quench into T = 0.5e. (c) An amorphous cluster obtained 
from the quench into T = O.Ole. 



(a) 



(b) 





Fig. 3. (Color online) The time evolution of an impact process of a LJ cluster of 300 atoms on the 
substrate, where the incident speed is Vimp = 2.0-^/ e/m. See the text in details. 



§3. Results 

3.1. Time evolution of impact processes 

Figures [3] and [H display the time evolutions of the impact process of the LJ 
cluster of 300 atoms on the crystalline surface. Figures [3] (a)-(d) represent the case 
of Mmp = 2.0y^ e/m, and Figs. [4] (a)-(d) are the case of ^imp = 4.0y^ e/m. 

The incident cluster moves toward the substrate with its translational speed 
Vimp ( Figs. El (a) and|l(a) ), and hits the substrate ( Figs. [3] (b) and[l(b) ). After 
the hitting, the cluster is only deformed to be a hemi-sphere ( Fig. [3] (c) ) for the 
small incident speed. If the incident speed is, however, larger than a critical value, 
the deposited cluster is split into many pieces ( Fig. [11(c) ). After time goes on, the 
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deposited cluster is adsorbed on the substrate and settles into the final configuration 
( Figs.[3](d) and [31(d) ). 

We observe that the impact process and the final configuration depend strongly 
on the incident speed Mmp- In the case of V^mp < 1.7y^ e/m, no atoms can escape 
from the cluster during the impact. By contrast, some atoms evaporate during 
the impact process for Mmp ^ l.Ty^e/m. If the incident speed is relatively small, 
the final configuration is a hemi-sphere on the substrate, as in the case of a partial 
wetting of a liquid droplet on a dry surface. The deformation is larger as the incident 
speed increases. Above Mmp = S.Sy^eTwi, the deposited cluster is completely split 
into fragments and the absorbed atoms on the substrate form a monolayer coverage. 
Above Vimp = 4.5ye/m, the deposited cluster is burst into fragments, and the 
absorbed coverage is no longer characterized by one cluster. 

At the moment of the impact, the temperature of the deposited cluster increases 
because the initial kinetic energy is transformed into internal motion. Then, the 
temperature decreases due to the heat conduction into the bulk region of the mate- 
rial through the contact area.'^''^'S3' The configuration of the deposited cluster 
is changed into an energetically favorite position during this cooling process. Fur- 
thermore, the atomic structure of the deposited cluster is adjusted to the substrate. 
Finally, the configuration is frozen because of the loss of temperature. 

3.2. The incident speed dependencies of the scaled cluster size and the cluster ad- 
sorption parameter 

In our simulation, the main cluster is detected by using the clustering algo- 
rithm.® Following the Allen and Tildesley, we adopt the critical atom separation 
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Fig. 5. (a) A plot of the scaled cluster size ^ and (b) a plot of the cluster adsorption parameter rj 
as the functions of the incident speed for the cluster of 300 atoms. 



Tc = 1.6(7. After the cluster settles into the final configuration, we represent N^is 
as the number of atoms in the cluster. With the aid of the number of atoms in the 
cluster before the impact N , we introduce the scaled cluster size : 



Nds 
N 



(3-1) 



If ^ = 1, no atoms can escape from the cluster after the impact. On the other hand, 
if < 1, evaporation of some atoms exists during the cluster impact. 

Let us define an absorbed atom in the cluster if an atom at r in the cluster 
satisfies the relation |r — r^l < r^, where is the position of its nearest neighbor 
constituent of the substrate. Using the number of these adsorbed atoms Nadh-, we 
can introduce the cluster adsorption parameter : 



•n 



Ncls 



(3-2) 



If ?7 < 1, the cluster is regarded as a multilayered adsorption. However, if = 1, 
the deposited cluster is perfectly spread on the substrate, and it is a monolayered 
adsorption. 

Figures [5] (a) and (b) plot the incident speed dependences of ^ and rj for the 
cluster of 300 atoms. We find that ^ equals to 1 below Vimp = l.Ty^ e/m, but it 
decreases above Mmp = l-Ty^ e/m. On the other hand, rj increases with the incident 
velocity below V^mp = 3.3y^ e/m, but it is satisfied to r/ ~ 1 above V^mp = S.Sy^TJm. 
Figure E] (a) plots several results on ^ for TV = 255, 300, 500 and 862. While Fig. E] 
(b) is r/ - ?7o for N = 32,108,255,300,500 and 862, where % is rj at Mmp = 0. It 
seems that — r/o is independent of the size of clusters, while exhibits weak size 
dependence. 

How can we understand the behaviors in Figs. [5] and [nf During the impact, the 
temperature in the cluster increases because the kinetic energy is transformed into 
internal motion.l^ We assume that the energy fiux to the substrate ^buik during 
the impact is small and the temperature becomes niaximuni Tijiax when the speed 
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Fig. 6. (Color online) (a) Plot of the scaled cluster size ^, and (b) plot of the scaled cluster 
adsorption parameter rj as the functions of the incident speed. 



of the center of mass of the cluster becomes zero. Thus, the energy conservation law 
can be written as 

^mNVlp + ^NTo - ^NT^ax + AS, (3-3) 

where Tq is the temperature of the cluster before the impact, and AS is the change 
of the surface energy. With the introduction of the surface tension 7, the height of 
the deposited cluster h, the contact radius of the deposited cluster R and the ratio 
(j) = h/R, AS is given by 

AS = -f{2TT(/)R^ - AttRI), (3-4) 

where Rq is the mean radius of the cluster before the impact. Introducing the mean 
area fraction of the contact area padh = ^adh/T^R^ ■, AS can be rewritten as 

AS = ^('^^-ATTRi\ . (3-5) 

V Padh J 



^ . 2AS , , 

To + y^^p-^, (3-6) 



From Eq. (j3-3p . T^ax satisfies 



Because the binding energy per atom in the cluster is roughly equal to e, the evapo- 
ration takes place at T^ax — £• In our simulation, Tq is much smaller than Tmax, and 
the last term on the right hand side of Eq. (I3-6P is negligible for large A^. Thus, the 



evaporation is believed to take place near Vimp — \/3e/m. In Figs. [5] (a) and [6] (a), 
the scaled cluster size becomes ^ < 1 above Vimp = 1.7ye/m, which is consistent 
with the above estimation. For the clusters with 32 and 108 atoms, decreases faster 
than the other cases. In such cases, we cannot ignore the last term on the right hand 
side of Eq. (l3^ . 

During the impact, an evaporated atom carries away the volume energy uy which 
is the potential energy per atom and the kinetic energy ^Tmax from the cluster. We 
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assume that the internal energy of the deposited cluster decreases because of the 
energy flux to the bulk of the substrate. Therefore, after the cluster settles into the 
final configuration, the energy conservation law can be written as 



\mNVl,^ + ^iVTo = AS + <Ptuik + (1 - ONE, (3-7) 



where E = ^T^ax — uy is the energy carried away by an evaporated atom. Here, 
(1 — ^)A^ represents the number of evaporated atoms. If the incident kinetic energy 
is not large, the number of atoms in the cluster is approximately preserved during 
the impact. Therefore, it is reasonable that the scaled cluster size satisfies ^ ~ 1. If 
we assume </> — 1, the energy conservation law Eq. (13-71) can be simplified as 



mpadh.r2 , Padh f 4:TT-fRo ^bulk , 3 \ 



where we have used Eqs. (|3-2p and (|3-5p . 

We use padh = 0.91, because the adsorbed atoms should match the hexagonal 
lattice on the substrate. The mean radius Rq of a cluster consisting atoms satisfies 
Rq = rQN^/^ where we use rg = 0.68cr as a fitting parameter. The solid line in 
Fig. [5] (b) is the theoretical prediction ( Eq. (j3-8p ), where the surface tension 7 ~ 
3.5e/cJ^ and the energy flux per atom ipbulk = ^bulk/N ~ 1.7e are other two fitting 
parameters. 

The second term on the right hand side of Eq. (|3-8p can be written as 

m{N) = ^ ^ATrrhN-'/' - ^,^1, + ^Tq) . (3-9) 

It is interesting that 77 — r]Q{N) is independent of the cluster size. Figure [6] (b) 
plots our numerical results r] — r]Q{N) for < 1, which support the validity of the 
theoretical prediction ( the solid line ). 

3.3. Transition from partial wetting to perfect wetting of the deposited cluster 

Let us introduce the radius of the equimolar dividing surface ( Gibbs Surface 



R^ = / -^--^rMr (3-10) 

Padh Jo dr 

as the contact radius of a deposited cluster, where p{r) is the area fraction of the 
contact area with radial distance from the center of mass of the adsorbed atoms in the 
cluster. We also define the cluster height h as z^ax — zq, where Zmax is the maximum 
vertical position in the atoms in the cluster, and zq is the minimum vertical position. 
Assuming a meniscus shape to the deposited cluster, we geometrically calculate the 
contact angle 6. 

Figure [7] displays h, R and 6 for the deposited cluster consisting of 300 atoms 
as the functions of the incident speed. We observe that the cluster height h de- 
creases and the contact radius R increases as the incident speed increases. Above 
^mp = 3.3y^ e/m, the height and the contact angle becomes zero, which implies 
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Yimpjmh Vmpym/8 

Fig. 7. (a) Plots of the height h ( filled circle ) and the contact radius 7? ( open square ) and (b) a 
plot of the contact angle 6 as the functions of the incident speed for the LJ cluster of 300 atoms. 



that the deposited cluster becomes a monolayer film and is perfectly wetting on the 
substrate. In this regime, the monolayer film is spread further and its boundary is 
partially chipped. Therefore the contact radius decreases. The clusters consisting 
of 108, 255, 500 and 864 atoms are also perfectly wetting on the substrate at critical 
incident velocities. However, it should be noted that the transition from multilayer 
film to monolayer film is only the morphological change of the deposited cluster. In 
the case of the cluster consisting of 32 atoms, the number of adsorbed atoms is too 
few to define the wetting parameters h, R and 6. 

3.4. Morphology of the final configuration of the adsorbed atoms in the cluster 
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Fig. 8. Configurations of the adsorbed atoms in a deposited cluster of 300 atoms for each incident 
speed. Vimp equals (a) 0.5\/e/m, (b) 1.0a/ e/m, (c) 1.5a/ e/m, (d) 2.0\/ e/m, (e) 2.5a/ e/m, (f) 
3.0-^e/m, (g) 3.5^ ^/m, and (h) 4.0i/e/m, respectively. 



The boundary shape of the contact area depends strongly on the incident speed. 
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In order to investigate the morphology of the boundary shape, we define the radial 
distance of the boundary r = f{ip). Here, r and ip are the usual radial and azimuthal 
coordinates. We take the origin to the center of mass of the adsorbed atoms in the 
cluster. Moreover, we define a dimensionless variable g{'ip) for the boundarjSU a,s 

fW - R 



R 



(3-11) 

gini/)^ with the integer n = 



We also use its Fourier representations g{ilj) 
0,±1,±2,---. 

Figures [8] (a)-(h) are the variety of the final horizontal configurations of the 
adsorbed atoms in the cluster consisting of 300 atoms. We find that there are three 
phases in the boundary shape. Below Vimp = 1.5ye/m, the boundary shape is grainy 
( Figs. [8] (a), (b) ), and l^^p has some peaks at higher modes. This may be caused by 
the small number of adsorbed atoms. In the case of 1.5ye/m < V^^p < 3.0ye/m, 
the deposited cluster is uniformly spread on the substrate ( Figs. [8] (c), (d), (e), 

(f) ). Thus, peaks of |(7nP vanish and their boundaries can be fitted by circles. In 
the case of 3.0y^e/m < T^mp < 4.0ye/m, the deposited cluster becomes a thin film 
or a monolayer film. In this regime, the boundary is partially chipped ( Figs. [8] 

(g) , (h) ), where |g'nP has intense peaks at some modes in our simulation. Above 
Vimp = 4.0ye/m, the deposited cluster is burst into fragments and the number of 
adsorbed atoms is too few to define R. Thus, we can not define g{il^) in this regime. 

The thermal fiuctuation of a circular geometry step is estimated as (IffnP) = 
T /2Ti(iR'n? from the equipartition of energy among the g-n^^^^ Here, (3 is the step 
edge stiffness. In our simulation, the thermal fiuctuation is estimated as (Is'nP) ~ 
0.1, while the \gn\'^ has peaks ranging from 1000 to 3000 above T^mp = 'i-^ye/m. 
Therefore this intensive peaks refiect on the growth of some unstable modes of gn 
during the spread of the deposited cluster on the substrate. 

3.5. Wettability between different Lennard-, J ones atoms 
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Fig. 9. (a) Plots of the cluster adsorption parameter ?y, and (b) plots of the contact angle as the 
functions of the incident speed for both A-C case ( open square ) and A- A case ( filled circle ). 



To investigate the influence of the interaction energy, we also perform the molec- 
ular dynamics simulation in which the LJ parameters of A and C are used for the 
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atoms in the cluster and the atoms of the substrate, respectively. Henceforth, we call 
this situation A-C case. The Lorentz-Berthelot rule Eq. (j2-2p is adopted to calculate 
the cross LJ parameters eAc and (Tac- Other simulation settings are the same as §[2j 
Figure [9] displays the cluster adsorption parameter -q and the contact angle 6 as 
the functions of the incident speed for the LJ cluster of 300 atoms. We also plot 
the data in the case of the LJ parameters of A-A case. We observe that the cluster 
adsorption parameter remains low value even when the impact speed is high ( 
Fig. [9] (a) ). We also stress that any A cluster does not become a monolayer film 
in which 9 becomes zero during the deposition onto a C surface. This is because 
the wetted state of the argon cluster is unfavorite on the crystalline carbon surface, 
which is resulted from eAc < £• Thus, it is clear that not only the incident velocity, 
but also the choice of composites is important to determine the final configuration 
of the deposited cluster. 

§4. Discussion and Conclusion 

In this paper, the incident kinetic energy per atom in the cluster is less than 2 eV. 
In this case, the damage of the substrate due to the impact of a cluster can be ignored 
and we considered a single layer substrate. However, the influence of the interaction 
between the deposited cluster and the bulk of the substrate should be important. 
In general, the adsorption state is strongly influenced by the surface temperature,'^ 
but the substrate was assumed to be at T = before the cluster impact in our 
simulation. Therefore, the influence of the surface temperature is also important for 
future study. Moreover, we performed simulation only for one deposition event of 
the cluster at each incident speed. Thus we should take ensemble average of impact 
processes for future study. 

In conclusion, we found that deposited LJ clusters consisting 32, 108, 255, 300, 
500 and 862 atoms exhibit a transition from multilayered adsorption to monolayered 
adsorption at the critical incident speed, Vimp = 3.3^/eJm. From our simulation, 
we clarified that the deposited clusters are perfectly wetting on the substrate above 
the critical incident speed. Employing the energy conservation law, we estimated 
the critical value of the incident speed at which the evaporation begin to occur 
during the impact. The estimated critical value, l^mp = l.Ty^e/m, exhibits good 
agreement with our simulation results of By using the energy conservation law, 
we also found that the scaled cluster adsorption parameter is independent of the 
cluster size and is proportional to V-^^. These results exhibit good agreement with 
our simulation results. We performed the Fourier analysis of g{ip) and found that 
some modes becomes unstable for 3.0y^e/m < Vi^p < A.Oy^e/m. We also performed 
the molecular dynamics simulation of A-C case and we found that the A cluster does 
not become a monolayer film on the C surface. Thus, we concluded that not only 
the incident speed, but also the strength of the interaction between the cluster and 
the substrate is important to form a monolayer film on a substrate. 
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